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Abstract 

We present results for the top quark pair cross section at the Tevatron and the LHC. 
We use the resummed double differential cross section, employing the fully kinematics- 
dependent soft anomalous dimension matrices, to calculate the soft-gluon contributions 
at next-to-next-to- leading order (NNLO). We improve and update our previous estimates 
by refining our methods, including further subleading terms, and employing the most 
recent parton distribution function sets. The NNLO soft corrections significantly enhance 
the NLO cross section while considerably reducing the scale dependence. We provide a 
detailed discussion of all theoretical uncertainties in our calculation, including kinematics, 
scale, and parton distributions uncertainties and clarify the differences between our work 
and other approaches in the literature. 



1 Introduction 



The top quark holds a special place in the Standard Model of particle physics as the heaviest 
elementary particle. Since the discovery of the top quark via tt production in proton-antiproton 
collisions at the Tevatron [T] its mass [2] and production cross section [3lll] have been determined 
with increasing accuracy. There is also now evidence for single-top production at the Tevatron 
[5]. At the LHC, both the ti and single top production cross sections will be two orders of 
magnitude higher than at the Tevatron. For recent reviews of top quark physics in hadron 
colliders see pj. 

The top quark cross section receives large corrections from soft-gluon contributions near 
threshold which can be formally resummed. The resummation at next-to-leading logarithmic 
(NLL) accuracy for hard scattering cross sections and, in particular, top quark pair production 
was performed in Ref. i?]. (For recent results on single top production, see Ref. |8j). To achieve 
such accuracy, it is necessary to derive the soft anomalous dimension matrix, which controls 
noncollinear soft-gluon emission, to one loop. At NLL the color structure of the hard scattering 
enters in a non-trivial way and each partonic process has to be treated separately. The soft 
anomalous dimension matrix is dependent on all the kinematical variables. Thus this is a fully 
differential calculation which can be applied to total cross sections as well as to differential 
cross sections, such as transverse momentum and rapidity distributions. 

Later, another formalism [5] was proposed for the total cross section only. This calcu- 
lationally simpler approach does not, however, involve the exact differential kinematics and 
instead makes the approximation that the NNLO and NLO rapidity dependence is the same. 
Hence, numerical deviations from the exact kinematics-sensitive result can appear. (For a de- 
tailed discussion and a numerical comparison in the context of direct photon production, see 
Ref. HOj). 

In Refs. O [9], the resummation is performed in moment space. Since the expression for 
the resummed cross section diverges at the Landau pole, a prescription is needed to define the 
physical resummed cross section when inverting from moment to momentum space. Alterna- 
tively, to avoid prescription dependence, the resummed cross section can be expanded to NNLO 
or higher orders. 

The formalism of Ref. [7| was used in detailed phenomenological studies [TT| |T2] where 
NNLO expansions were provided at NNLL accuracy, after matching with the complete NLO 
cross section. Results were provided in both single-particle-inclusive (IPI) and pair-invariant- 
mass (PIM) kinematics. The kinematics ambiguity was found to be an important source of 
uncertainty. In IPI kinematics the soft logarithms are of the form [In^ (34/711?) / s 4]^ with m the 
top quark mass and S4 the sum of the Mandelstam invariants, S4 = s + ti +ui. Near threshold, 
S4 — * 0. The soft-gluon corrections to the double differential cross section, (fa/ {dtidui), were 
calculated. In PIM kinematics, the soft logarithms are of the form [ln'^(l — z)/{l — z)]^ with 
z = M'^/s, where is the ti pair mass squared. Near threshold, z ^ 1. The soft gluon 
corrections to the double differential cross section, (P a /{dM'^d cos 6), where 6 is the scattering 
angle in the partonic center-of-mass frame, were calculated. The cross section in PIM kinematics 
was found to be smaller than the IPI result. The difference, an uncertainty due to uncalculated 
terms, was found to be larger than the scale variation. In Ref. p^], results were also given for 
the exact scale variation at NNLO. The magnitude of this variation was also found to depend 
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on the kinematics. 

The formahsm of Ref. [9] was also used to derive NLL resummed numerical results (recently 
updated [13]) for top pair production at the Tevatron and LHC The corrections beyond NNLO 
are negligible [S] (also shown to be small in Ref. [E]); hence the resummed cross section is 
numerically very similar to the NNLO expansion at the given logarithmic accuracy. However, 
a minimal prescription was used to define the resummed cross section in Ref. [9] and, as shown 
in Ref. [11], the differences in the prescription formalism, as well as the treatment of the 
kinematics, are much bigger than higher order terms at NNNLO and beyond. Hence the 
results of Ref. [9] are quite different from Refs. [TTl [12] , both theoretically and numerically. 

The approach in Refs. [TT], [12] was later improved by adding NNNLL terms at NNLO 
|15j . Although the complete NNNLL terms require calculation of the two- loop soft anomalous 
dimension matrix, it was clearly demonstrated [15] that the contribution of this matrix at 
two loops is expected to be negligible. Thus it is possible to obtain an effective NNNLL 
calculation by including all other terms. The ( terms arising from the inversion to momentum 
space (including some ( virtual terms) are dominant, as shown by expressing the partonic 
cross sections in terms of scaling functions that depend on the variable rj = s/(4m^) — 1 and 
comparing the IPI and PIM scaling functions over a large range of t]. Since a complete NNLO 
calculation should be independent of the kinematics, the difference between the IPI and PIM 
results as a function of rj is an indication of the unknown terms. Away from threshold, hard 
gluon terms contribute. Since their form is also kinematics dependent it is inevitable that, as 
one moves away from the threshold region, the IPI and PIM results diverge. However, near 
threshold the soft gluons dominate and thus a complete calculation of the soft terms should 
produce agreement between the IPI and PIM scaling functions. At NNLL, the IPI and PIM 
functions diverge already at threshold [12j, indicating that the NNNLL terms are non-negligible. 
However, when the NNNLL terms were added [TJ], this discrepancy disappeared, as one would 
expect when all NNLO soft terms are included. Thus the contribution of the unknown two- 
loop soft-anomalous-dimension terms that were left out is negligible and we obtain an excellent 
approximation to the complete NNNLL terms and an effective NNNLL calculation, denoted 
NNLO-NNNLL+C in Ref. [I5]. 

Recently the two-loop soft anomalous dimension for massless quark scattering was completed 
|16j . (Work is in progress for heavy quark production [I?].) It was shown [TB] that the two-loop 
soft anomalous dimension is simply the one-loop result multiplied by half the two-loop quantity 
K [18]. Assuming that this relation also holds for heavy quarks, consistent with the two-loop 
results in Ref. [19], the contribution of this additional two- loop term to the total cross section 
is less than 1 per mille at both the Tevatron and LHC energies. It is thus insignificant relative 
to the size of other terms and sources of uncertainty, as expected [15], verifying the robustness 
of the calculation in Ref. [T5J. 

A very recent paper [19] uses the general approach of Ref. [9J, extending the results of 
Ref. [9j by adding the NNLL terms in the resummed expression. A NNLO expansion in powers 

of In P, where /? = y'l — Am? / is also presented in Ref. [19j . An additional two-loop term is 
also included, a rough analog of the two-loop soft anomalous dimension term in the formalism 
of Refs. [H [ini [IZ]. This two- loop term is again given by the one-loop result multiplied by 
the two-loop quantity K [IB], analogous to the result of Ref. [IB]. We have investigated the 
contribution of this two-loop term within the approach of Ref. [I9] and find it to be numerically 
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negligible, on the order of a few per mille, consistent with the study mentioned previously and 
again verifying that this two-loop contribution is numerically insignificant. Since [19] uses the 
formalism and approximation of Ref . [9j , their results differ from ours for the reasons explained 
earlier. We use exact kinematics in a double differential cross section and define partonic 
threshold through the quantities S4 or z, depending on the kinematics choice, while Ref. [T^ 
defines threshold at the total cross section level only in terms of (3. 

We also note that the authors of Ref. [19] use the exact scale dependence of the NNLO 
cross section, finding a smaller scale dependence than in Refs. [9l [13] . The exact NNLO scale 
dependence was first calculated in Ref. [12] and shown to crucially depend on the kinematics 
choice (IPI or PIM). Hence the scale variation in Refs. [9], [19] cannot be directly compared 
with that of Ref. [12j or the present paper. 

We also find it more consistent to use the same level of accuracy for the scale-dependent 
terms as for the other terms in the calculation, also chosen for the final results in Refs. [T2| [T5]. 
We obtain a smaller scale dependence at the Tevatron than Ref. [12] but a larger one at the 
LHC In our approach, the kinematics ambiguity is bigger than the scale variation at the 
Tevatron but smaller at the LHC. The results in Refs. [5], |T3l [T^] do not have a kinematics 
uncertainty because their approximation is insensitive to the kinematics choice. Therefore the 
scale dependence in those approaches, in particular the small scale uncertainty at the LHC of 
Ref. [19], is not necessarily indicative of the true theoretical uncertainty in the cross section. 

In Ref. [19] subleading Coulomb terms, calculated in [20], were included in the numerical 
results. In Ref. [15] some, but not all, of these terms were included. We find that these 
additional contributions are completely negligible at the Tevatron, and make a very small 
contribution, included in our new results, at the LHC. 

In this paper, we present detailed results for the top quark cross section at the Tevatron 
and the LHC. We primarily use the theoretical approach of Ref. [TJ] with some changes and 
refinements, described in the text, along with the newest available parton distribution functions 
(PDFs). We also provide a detailed study of theoretical uncertainties including kinematics, 
scale, and PDF uncertainties, as well as a discussion of other sources. In Section 2 we provide 
results for the ti production cross section in pp collisions at the Tevatron at =1.96 TeV. 
In Section 3 we give the ti production cross section in pp collisions at the LHC at = 14 
TeV as well as a prediction for 10 TeV. 

2 The top quark cross section at the Tevatron 

We begin with ti production at the Tevatron at = 1.96 TeV. The leading-order partonic 
processes are qq — ^ ti and gg — > ti. In pp collisions at the Tevatron, the qq channel is dominant. 
In addition to corrections to the LO processes, at NLO there are small contributions from two 
additional processes, qg and qg. 

We first calculate the NLO cross section [2T] , including all channels. We then add the NNLO 
soft-gluon corrections in the qq and gg channels to the NLO result. We calculate the soft-gluon 
corrections in both IPI and PIM kinematics. We find that the behavior of the qq — *■ ti and 
gg — >■ ti contributions is quite different. The qq channel, with only one diagram at LO, is 
well behaved in both kinematics. The gg channel is, however, better treated in IPI kinematics 
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because the three LO gg tt diagrams favor IPI kinematics. In addition, in PIM kinematics 
the one-loop expansion of the resummed gg cross section is very different from the exact NLO 
resuh while the IPI expansion of the gg contribution is an excellent approximation to the exact 
NLO result at both the Tevatron and the LHC Therefore, for our best prediction of the NNLO 
corrections, we take the average of the IPI and PIM qq soft-gluon corrections and add it to the 
IPI gg soft-gluon result. This is a refinement of our previous method [15] where we averaged 
the IPI and PIM results from both channels. Our new approach results in a slightly larger total 
cross section with a somewhat reduced kinematics uncertainty at the Tevatron than in Ref. [15]. 
As discussed in the Introduction, the effect of the two-loop soft anomalous dimension matrix 
and further subleading Coulomb terms is negligible with no change on the results presented 
here. 

We present two tables with NLO and approximate NNLO cross sections using the MRST 
2006 NNLO [22] and CTEQ6.6M [23] parton densities. While the NNLO approximate cross 
section is of the same logarithmic accuracy as the NNLO-NNNLL-|-(^ cross section of Ref. [15] . 
they differ slightly because of the changes and refinements discussed above. All the results are 
in the MS scheme. 



MRST 2006 NNLO 


Mass 
(GeV) 


(j(NLO ± scale 
± PDF) (pb) 


cr(NNLO approx ± kinematics 
± scale ± PDF) (pb) 


165 


Q 90 +0.59 
-'■^'J -1.09 -0.23 


9.80 ± 0.38 


166 


Qo +U.5V +U.2V 
-1.06 -0.22 


9.48 ± 0.37 


167 


«p; +0.55 +0.26 
-1.02 -0.21 


9.17 ±0.36int 


168 


07 +0.53 +0.25 
O-'J' -0.99 -0.20 


8.88 ± 0.35 1-1 l-i 


169 


1 1 +0.51 +0.24 
-0.96 -0.19 


8.60 ± 0.34 


170 


7 OCT +0.50 +0.23 
' -o^ -0.93 -0.19 


O.OZ in U.OO _Q 29 _o.20 


171 


7 fin +U.23 
'•OU -0.90 -0.18 


8.06 ± 0.32 i^l 


172 


7 oe +0.46 +0.22 
'•<J0 -0.87 -0.18 


7.80 ± 0.31 tl^i tin 


173 


7 1 Q +0.45 +0.21 
'•-LO -0.84 -0.17 


7.56 ± 0.30 t^:^ t^-^ 


174 


a Qi +0.44 +0.20 
"•yj- -0.82 -0.17 


7.32 ± 0.29 tVA 


175 


R 7n +0.42 +0.19 
^■'^ -0.79 -0.16 


7.09 ± 0.28 


176 


+0.41 +0.19 
0-4y -0.77 -0.15 


6.87 ± 0.27 t',f, 


177 


^ 9Q +0.39 +0.18 
""^y -0.74 -0.15 


6.66 ± 0.26 


178 


^ in +0.38 +0.18 
"•-LU _o.72 -0.14 


6.46 ± 0.26 


179 


f- Qi +0.37 +0.17 
o-yj- -0.70 -0.14 


6.26 ± 0.25 ^l^i 


180 


f- 70 +0.36 +0.17 
o- ' -0.68 -0.13 


6.07 ± 0.24 



Table 1: The tt production cross section in pp collisions at the Tevatron with \fS = 1.96 TeV 
using the MRST 2006 NNLO PDFs. The exact NLO results are shown with the scale and 
PDF uncertainties while the approximate NNLO results include kinematics, scale, and PDF 
uncertainties. 

Table [T] provides the pp —>■ tt cross section for 165 < m < 180 GeV, in 1 GeV increments 
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calculated with the MRST 2006 NNLO PDFs. We give both the exact NLO and the approxi- 
mate NNLO cross sections. The central values are calculated with the factorization scale, hf, 
and the renormalization scale, fiR, set equal to the top quark mass, fip = fin = fi = m, using 
the central MRST 2006 NNLO PDF. 

In addition to the central value of the NLO cross section, we also provide uncertainties due 
to the scale and PDF variations. The scale is varied over m/2 < fi < 2m. Varying the scale by 
a factor of two around yU = m is a standard but arbitrary way to estimate uncertainties from 
higher-order terms. The + (— ) indicates the difference between the calculation with /i = m/2 
[fx = 2m) and the central value with fi = m. The PDF uncertainty, calculated using the 30 
different MRST 2006 NNLO eigensets, is relatively large, reflecting the uncertainty in the large 
X region of the PDFs. 

The central value of the NNLO approximate cross section is followed by the kinematics 
uncertainty, the scale variation and the PDF uncertainty. The NNLO scale and PDF uncer- 
tainties are obtained the same way as for the NLO cross section. The kinematics uncertainty 
resides in the treatment of the qq channel because the gg channel is only calculated in IPI 
kinematics. The central value is obtained from the average of the IPI and PIM qq calculations. 
The -|- kinematics uncertainty is the found by taking the IPI qq result alone while the — results 
uses the PIM qq calculation. The kinematics uncertainty is symmetric becuase the central qq 
contribution is the average of the IPI and PIM results. 

At NLO, the scale variation is significant. When the NNLO corrections are added, the scale 
dependence on the scale decreases dramatically. However, the kinematics uncertainty is larger 
than the scale variation. The PDF uncertainty is also significant at NLO and NNLO, of the 
same order as the kinematics dependence at NNLO. 

Table [2] provides the ti cross section for 165 < m < 180 GeV, in 1 GeV increments, using 
the CTEQ6.6M NLO PDFs. As in Table 1, we list both the exact NLO and our approximate 
NNLO cross sections together with all uncertainties. The central values, again shown with 
fJ'F = fJ'R = fJ' = "tTT-, employ the central CTEQ6.6M PDFs. 

The scale variation is calculated as described above for the MRST 2006 NNLO PDFs. In 
the case of the NNLO approximate cross section, the NNLO results at both ends of the scale 
range, /i = m/2 and 2m, are lower than with /z = m, indicated by the double minus signs on 
the scale uncertainty. The PDF uncertainty is calculated using the 44 different CTEQ6.6M 
eigensets. 

The two sets of results are quite different. The cross sections calculated with CTEQ6.6M are 
smaller than those with the central MRST 2006 NNLO set but have larger PDF uncertainties. 
Indeed, the CTEQ6.6M NNLO approximate cross sections are quite similar to the NLO cross 
section calculated with the MRST 2006 NNLO PDFs. Part of the difference can be attributed 
to the fact that the MRST 2006 NNLO sets are of the same order as the NNLO approximate 
calculation while the CTEQ6.6M sets are NLO. Furthermore, the CTEQ6.6M large-x gluon 
distribution is smaller, reducing the relative gg contribution. 

The best way to combine the uncertainties is not obvious. The most conservative approach 
would be to add them linearly. However the kinematics and scale uncertainties both reflect 
the neglect of uknown terms. Thus a linear combination of the uncertainties likely provides an 
overestimate of the overall uncertainty and we instead prefer to add them in quadrature. 

We present the NNLO approximate cross section for the current most likely value of the top 
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CTEQ6.6M 


Mass 
(GeV) 


a(NLO ± scale 
± PDF) (pb) 


cr(NNLO approx ± kinematics 
± scale ± PDF) (pb) 


165 


o 7/1 +U.46 +U.58 
o- ' ^ -0.96 -0.45 


Q 4- n 'iV ^0.03 +U.6i 
J.ZO 31 U.O / _o.25 -0.48 


166 


o /I 7 +0.44 +0.56 
o-^' -0.93 -0.43 


8.93 ± 0.36 -Jof, tnt 


167 


o on +'J-43 +0.53 
-0.90 -0.42 


8.65 ± 0.35 _o;23 


168 


7 c\A 

' ■■J^ -0.87 -0.41 


8.38 ± 0.34 Hit 


169 


7 70 +0.40 +0.50 
'•'U -0.84 -0.39 


010 -1- n 'i'i ^0.03 +U.53 

O.iZ HZ U.OO _o.22 -0.41 


170 


7 /ic +0.39 +0.48 
' -^"J -0.82 -0.38 


7.87 ± 0.32 till, 


171 


7 OQ +0.38 +0.47 
' -^"^ -0.79 -0.36 


7.62 ±0.311:^^^:^^ 


172 


7 ni +0.37 +0.45 
'•UJ- -0.77 -0.35 


7.39 ± 0.30 


173 


7Q +0.35 +0.43 
' -0.74 -0.34 


7.16 ±0.29 1^^+^^^ 


174 


c CO +0.34 +0.42 
O-Oo _o.72 -0.33 


6.94 ± 0.28 I^^ t^M 


175 


^ oo +0.33 +0.4i 
O-'JO -0.70 -0.31 


6.73 ± 0.27 Z'.Ts 


176 


a 1 Q +0.32 +0.39 
'^■^'J -0.68 -0.30 


0.00 =■= U./; ^ _o.i8 -0.32 


177 


nn +"-^1 +"-y« 

"■UU _o.66 -0.29 


6.33 ± 0.26 -J]-^? +1]-^'/ 


178 


cr 09 +0.30 +0.37 
^•o^ -0.64 -0.28 


6.14 ±0.25 Z^^H^ 


179 


p; «p; +0.29 +0.35 
-0.62 -0.27 


5.95 ± 0.24 Z',f, i^:^^ 


180 


c- /,o +0.28 +0.34 
'-'■^0 -0.60 -0.26 


K 77 _L n 9/1 -0.03 +0.36 
0.1 1 IC U.Z^ _Q _Q 27 



Table 2: The production cross section in pp collisions at the Tevatron with = 1.96 TeV 
using the CTEQ6.6M PDFs. The exact NLO results are shown with the scale and PDF uncer- 
tainties while the approximate NNLO results include kinematics, scale, and PDF uncertainties. 



quark mass, m = 172 GeV for both sets of PDFs. We first give the central result from Tables [T] 
and [2] with all the uncertainties shown separately and then add the uncertainties in quadrature 
for our final resuh. Using the MRST 2006 NNLO PDFs, we have 

a^''^,?"PP™"(1.96 TeV, m = 172 GeV, MRST) = 7.80 ± 0.31 +° °? t'oH pb = 7.80 pb , 

(2.1) 

while with CTEQ6.6M we find 
^NNLOapprox^^ gg rj.^y ^ ^ ^ ^ CTEQ) = 7.39 ±0.30 pb = 7.39 pb . (2.2) 

We have not included further theoretical ambiguities arising from the choice of equivalent 
analytical expressions near threshold or from damping factors [12j as well as from the virtual ( 
terms discussed in Ref. ^5]. Such ambiguities are partly accounted for in the kinematics and 
scale uncertainties shown here. 

Figure [1] shows the exact NLO and approximate NNLO top quark cross sections as a function 
of top quark mass at the Tevatron using the MRST 2006 NNLO (left) and CTEQ6.6M (right) 
PDFs. Three curves are given for each order: a central value with fi = m and the extremes of 
the calculated scale dependence with fi = m/2 and 2m. The region between the upper and lower 
scales represents the scale variation at each order. We see that the NNLO scale dependence is 
much diminished relative to NLO. In fact the NNLO curves with ^ = m/2 and fi = m are on 
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Figure 1: The exact NLO and approximate NNLO cross sections in pp collisions at 1.96 TeV 
using the MRST 2006 NNLO (left) and CTEQ6.6M (right) PDFs. 



top of each other as well as on top of the NLO curve with n = m/2. The kinematics and PDF 
uncertainties are not represented in the plots. 

In Fig. El we present the K factors at the Tevatron using the MRST 2006 NNLO (left) 
and CTEQ6.6M (right) PDFs with our central value, /i = m. The K factors are virtually 
independent of the PDF although the CTEQ6.6M K factors appear slightly smaller. They are 
also independent of the top quark mass. The ratios of the approximate NNLO cross sections 
to the exact LO and NLO cross sections are both given. The NNLO corrections enhance the 
NLO ti cross section by ~ 6% at /i = m. 

3 The top quark cross section at the LHC 

We now turn to ti production in pp collisions at the LHC. While our results are primarily 
shown for \fS = 14 TeV, we also provide predictions for the top quark cross section at the LHC 
start-up energy of 10 TeV. We note that at the LHC, the gg channel is dominant. 

Table [3] provides the top quark cross section for 165 < m < 180 GeV, in 1 GeV increments, 
in pp collisions at the LHC with = 14 TeV employing the MRST 2006 NNLO PDFs. We 
list both the exact NLO and approximate NNLO cross sections. The central values are given 
for /iiT = /iij = /i = m with the central MRST 2006 NNLO PDFs. 

The NLO cross section is shown with the uncertainties due to the scale variation and the 
choice of PDF eigenset. As before, we vary the scale between /i = m/2 and 2m and the PDF 
uncertainty is calculated using the 30 different MRST 2006 NNLO eigensets. 

The central value of the NNLO approximate cross section is accompanied by uncertainties 
due to the kinematics, the scale variation and the choice of PDF. Again, the central value of 
the approximate NNLO qq contribution if the average of the IPI and PIM kinematics choice. 
The + (— ) kinematics uncertainty is the difference between the top cross section with the qq 
contribution calculated in IPI (PIM) kinematics. 
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pp->tt Tevatron S =1.96TeV n=m MRST2006 pdf 



pp->tt Tevatron S =1.96TeV n=m CTEQ6.6 pdf 



1 


1 1 1 1 






NLO/LO 

NNLO approx / LO 

NNLO approx / NLO 


















1 1 1 1 1 1 1 1 1 1 1 1 1 1 



170 175 
m (GeV) 



1.6 



1.2 



165 



NLO/LO 

NNLO approx / LO 

NNLO approx / NLO 



_l I 1_ 



170 175 
m (GeV) 



_l I 1_ 



180 



Figure 2: The K factors at the Tevatron using the MRST 2006 NNLO (left) and the CTEQ6.6M 
(right) PDFs. 



At NLO the scale variation is large. When the NNLO corrections are added the dependence 
on the scale decreases significantly. Now the kinematics uncertainty is much smaller than that 
due to the scale variation. This is because we only use IPI kinematics for the gg channel, 
dominant for pp collisions at this energy. Thus the kinematics uncertainty is only due to the 
change in the qq calculation. The PDF uncertainty is smaller at the LHC since x is relatively 
small, in a range where the PDFs are better known. 

In TablelU we present the corresponding top cross sections with the CTEQ6.6M NLO PDFs. 
The CTEQ6.6M results are again smaller than those with the MRST 2006 NNLO sets. While 
the PDF uncertainty is smaller at the LHC, the CTEQ6.6 uncertainty is still larger than those 
with MRST 2006 NNLO. 

We now present our predicted NNLO approximate cross section for top production in pp 
collisions at = 14 TeV with m = 172 GeV. The results are again given both with the 
separate uncertainties, as in Tables [3] and HJ and with uncertainties added in quadrature. Using 
the MRST 2006 NNLO PDFs, we find 

%'!!;?''''"■"( 14 TeV, m = 172 GeV, MRST) = 968 ±4 pb = 968 pb , (3.1) 

while with the CTEQ6.6M PDFs, we obtain 

%'!I;?PP''°'(14 TeV, m = 172 GeV, CTEQ) = 919 ±4 1™ pb = 919 tH pb . (3.2) 

Figure [3] shows exact NLO and approximate NNLO top quark cross sections as a function of 
top quark mass with = 14 TeV using the MRST 2006 NNLO (left) and CTEQ6.6M (right) 
PDFs. At each order we show the central result with yU = m as well as the range of the scale 
uncertainty indicated by the upper {fi = m/2) and lower (/i = 2m) curves. The region between 
the upper and lower curves denotes the scale variation. While the NNLO scale dependence 
is reduced relative to NLO, the reduction is not as large as at the Tevatron. We note that 
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pp->tt at LHC S =14TeV MRST2006 pdf 



1400 r 



— 1 1 1 1 1 1 1 — 



NLO n=in 
NLO n=m/2, 2m 
NNLO approx [i=m 
NNLO approx |i=m/2, 2m 



3 1000 ■ 
D 



800 - 



pp->tt at LHC S =14TeV CTEQ6.6 pdf 



1400 



& 1000 



800 



— 1 — I — I — 



— 1 1 r- 



NLO |l=m 
NLO ^l=m/2, 2m 
NNLO approx [l=m 
NNLO approx n=m/2, 2m 



^65 170 175 180 *''?65 170 175 180 

m (GeV) m (GeV) 

Figure 3: The NLO and approximate NNLO top cross sections in 14 TeV pp collisions at the 
LHC using the MRST 2006 NNLO (left) and the CTEQ6.6M (right) PDFs. 



pp->tt LHC S =14 TeV n=m MRST2006 pdf 



pp->tt LHC S =14 TeV [i=m CTEQ6.6 pdf 
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NNLO approx / LO 
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NLO/LO 
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Figure 4: The K factors at the LHC using the MRST 2006 NNLO (left) and CTEQ6.6M (right) 
PDFs. 
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MRST 2006 NNLO 



Mass 
(GeV) 


a(NLO ± scale 
± PDF) (pb) 


cr(NNLO approx ± kinematics 
± scale ± PDF) (pb) 


165 


1089 l[i 


1173 ± 5 


166 


1059 +li t[i 


1141 ± 5 rti 1^1 


167 


1030 t[i 


1109 ± 5 tf^ 11^ 


168 


1003 t[i 


1080 ± 5 


169 


976 1- 


1050 ± 5 til 


170 


950 li- 1- 


1022 ± 5 tf, t'^ 


171 


Q9/I +ii4 

_no -12 


994 ± 5 11^ l^i 


172 


900 ^^^"^ 

JUU _io7 -12 


968 ± 4 


173 


07« +iU8 +11 
O'O -104 -11 


943 ± 4 


174 


Qco +1U5 +11 
-101 -11 


918 ±4^^^^!^ 


175 


001 +1U2 +11 
""JJ- -98 -11 


894 ±AV^ t^i 


176 


ono +yy +10 

8Uy _96 _io 


871 ± 4 


177 


7QO +yy +1U 

'^^ -93 -10 


848 ±AtZt{{ 


178 


7R8 +^4 +1U 
'Oo -91 -10 


826 ± 4 t{{ 


179 


748 

-89 -10 


805 ± 4 +ii 


180 


79Q +89 +9 
-86 -10 


785 ± 4 1^^ t^i 



Table 3: The ti production cross section in pp collisions at the LHC with = 14 TeV using the 
MRST 2006 NNLO PDFs. The exact NLO results are shown with scale and PDF uncertainties 
while the approximate NNLO results include kinematics, scale, and PDF uncertainties. 



the mass dependence at the LHC is smaller than at the Tevatron since we are further from 
production threshold here. 

In Fig. m we show the LHC K factors as a function of mass for our central (/i = m) cross 
sections calculated with the MRST 2006 NNLO (left) and CTEQ6.6M (right) PDFs. The LHC 
K factors are larger than those shown in Fig. O They are almost identical for the two sets and 
are virtually independent of mass. The approximate NNLO cross section is ~ 8% larger than 
the NLO cross section. 

Finally, we provide predictions for the initial LHC run at = 10 TeV with m = 172 GeV. 
Using the MRST 2006 NNLO PDFs, the exact NLO cross section is 414 ±52 ± 8 pb while 
the NNLO approximate cross section is 

V'!!;?''^^™"(10 TeV, m = 172 GeV, MRST) = 446 ±3 1| ± 9 pb = 446 1| pb . (3.3) 

With the CTEQ6.6M PDFs we find that the NLO cross section is 385 t% tJl pb and the NNLO 
approximate cross section is 

V'!l^i?"''''™'(10 TeV, m = 172 GeV, CTEQ) =415 ±2 tfi ± 20 pb = 415 pb . (3.4) 
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CTEQ6.6M 



Mass 
(GeV) 


cr(NLO ± scale 
± PDF) (pb) 


(t(NNLO approx ± kinematics 
± scale ± PDF) (pb) 


165 


lUOO _i2i _34 


1114 ±5-^ ^r, 


166 


1007 t^ii ±i 


1084 ± 5 ttl 


167 


979 lii^ 


1054 ± 5 tti li 


168 


y^-^ -112 -32 


1025 ± 5 ^^,1 If, 


169 


Q97 +110 +29 
-109 -31 


997 ± 4 mi 


170 


Qno +107 +29 

_IQQ _3Q 


970 ± 4 1- ^- 


171 


Q77 +105 +28 
o' ' -103 -29 


944 ± 4 


172 


or 1 +H)2 -2i 
Od l _]^00 -29 


919 ± 1 -11^ 


173 


oqi +99 +27 
OOi _29 


894 ± 4 +^5* 


174 


onn +t»ti +^ti 
oUy _95 _28 


870 ±4^^^!^^ 


175 


700 +94 +26 
'00 -92 -27 


847 ± 4 


176 


7/^7 +91 +25 
-90 -27 


825 ±41^^!^^ 


177 


7/1 7 +89 +25 

-88 -26 


803 ± 4 1^;; 


178 


707 +86 +25 

-85 -26 


782 ±41^11 ii 


179 


709 tit li 


762 ± 3 


180 


690 1^1 -25 


742 ± 3 



Table 4: The production cross section in pp collisions at the LHC with y/S — 14 TeV using 
the CTEQ6.6M PDFs. The exact NLO results arc shown with scale and PDF uncertainties 
while the approximate NNLO results include kinematics, scale, and PDF uncertainties. 

4 Conclusions 

We have studied top quark production at the Tevatron and the LHC. Our work is the only 
calculation that employs full kinematics in the double differential cross section beyond NLL 
using the soft anomalous dimension matrix. We presented detailed results for the exact NLO 
and approximate NNLO tt cross sections at the Tevatron and the LHC for a wide range of 
top quark masses using the MRST 2006 NNLO and the CTEQ6.6M PDFs. The approximate 
NNLO corrections include soft-gluon contributions which significantly enhance the cross section. 
We find that further two-loop soft gluon contributions are expected to be negligible. We also 
included subleading Coulomb contributions and found them negligible at the Tevatron and very 
small at the LHC. 

We provided detailed results for the theoretical uncertainties, including the kinematics am- 
biguity, scale variation, and PDF uncertainties. We found that the NNLO scale uncertainty is 
drastically reduced relative to NLO at the Tevatron where the kinematics uncertainty is larger. 
The PDF uncertainty is quite significant, especially for the CTEQ6.6M PDFs. At the LHC, the 
kinematics ambiguity is small. The NNLO scale variation is larger despite being significantly 
smaller than at NLO. The PDF uncertainty is smaller at the LHC than at the Tevatron. The 
results using the MRST 2006 NNLO and CTEQ6.6M PDFs are quite different fi-om each other 
at both Tevatron and LHC energies. 
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Ongoing work with two-loop soft anomalous dimensions in the eikonal approximation [T7j 
and recent analytical two-loop pieces of the NNLO corrections [23] promise further progress in 
the future. 



Acknowledgements 

The work of N.K. was supported by the National Science Foundation under Grant No. PHY 
0555372. The work of R.V. was performed under the auspices of the U.S. Department of 
Energy by Lawrence Livermore National Security, LLC, Lawrence Livermore National Labora- 
tory under Contract DE-AC52-07NA27344 and also supported in part by the National Science 
Foundation Grant NSF PHY-0555660. 



References 

[1] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 74, 2626 (1995) |hep-ex /9503002|; 
DO Collaboration, S. Abachi et al, Phys. Rev. Lett. 74, 2632 (1995) |hep^ x/9503003|. 

[2] Tevatron Electroweak Working Group for CDF and DO. l arXiv:0803.16 83 [hep-ex]. 

[3] CDF Collaboration, Phys. Rev. Lett. 96, 202002 (2006) |hep-ex/0603043l; Phys. Rev. D 
74, 072006 (2006) |hep-ex/0607035j; Phys. Rev. D 74, 072005 (2006) | hep-ex/0607095j ; 
Phys. Rev. D 76, 072009 (2007), larXiv:0706.3790 [hep-ex]. 

[4] DO Collaboration, Phys. Rev. D 74, 112004 (2006) |hep-ex/0611002]; Phys. Rev. D 76, 
072007 (2007) |hep-ex/0612040|; Phys. Rev. D 76, 092007 (2007), arXiv:0705.2788 [hep- 
ex]; Phys. Rev. D 76, 052006 (2007). larXi v:0706 .0458 [hep-ex]; arXiv:0803.2779 [hep-ex]. 

[5] DO Collaboration, V.M. Abazov et al, Phys. Rev. Lett. 98, 181802 (2007) 
|hep-ex/0612052|; arXiv:0803.0739 [hep-ex]; 
CDF Collaboration, Conf. Note 8964; Conf. Note 8968. 



[6] W. Wagner, Rept. Prog. Phys. 68, 2409 (2005) |hep-ph/0507207|; A. Quadt, Eur. Phys. 
J. C 48, 835 (2006); R. Kehoe, M. Narain, and A. Kumar, Int. J. Mod. Phys. A 23, 
353 r2008). [a^:0712.2733l [hep-ex]: T. Han. larXiv:0804.3178l [hep-ph]: W. Bernreuther, 
larXiv:0805. 13331 [hep-ph]. 

[7] N. Kid onakis and G. St erman, Phys. Lett. B 387, 867 (1996); Nucl. Phys. B505, 321 
(1997) |hep-ph/9705234j. 

[8] N. Kidonakis, Ph ys. R ev. D 74, 114012 (2006) (hi^h/0609287] ; Phys. Rev. D 75, 
071501(R) (2007) p ^ph/070108 0]. 

[9] R. Bonciani, S. Catani, M.L. Mangano, and P. Nason, Nucl. Phys. B529, 424 (1998) 
| hep-ph/980137 5]. 

[10] G. Sterman and W. Vogelsang, JHEP 02, 016 (2001) |hep-ph/0011289] . 



13 



[11] N. Kidonakis, Phys. Rev. D 64, 014009 (2001) |hep-ph/0010002] . 

[12] N. Kidonakis, E . Laenen, S. Moch, and R. Vogt, Phys. Rev. D 64, 114001 (2001) 
|hep-ph /0105041 ). 

[13] M. Cacciari, S. Frixione, M.L. Mangano, P. Nason, and G. Ridolfi, arXiv:0804.2800 [hep- 
ph]. 

[14] N. Kidonakis, Phys. Rev. D 73, 034001 (2006) |hep-ph/0509079] . 

[15] N. Kidonakis and R. Vogt, Phys. Rev. D 68, 114014 (2003) |hep-ph /0308222]. 

[16] S .M. A ybat, L.J. Dixon, and G. Sterman, P hys. Rev . Lett. 97, 072001 (2006) 
|hep-ph/0606254]; Phys. Rev. D74, 074004 (2006) |hep-ph/06 07309|. 

[17] N. Kidonakis and P. Stephens, in DIS 2008. larXiv: 0805. 1193 1 [hep-ph]. 

[18] J. Kodaira and L. Trentadue, Phys. Lett. 112B, 66 (1982). 

[19] S. Moch and P. Uwer, larXiv: 0804. 14761 [hep-ph]. 

[20] W. Bernreuther, R. Bonciani, T. Gehrmann, R. Heinesch, T. Leineweber, P. Mastroha, 
and E. Remiddi, Nucl. Phys. B706, 245 (2005) [hep-ph/0406046] . 

[21] P. Nason, S. Dawson, and R.K. Elhs, Nucl. Phys. B303, 607 (1988); 

W. Beenakker, H. Kuijf, W.L. van Neerven, and J. Smith, Phys. Rev. D 40, 54 (1989); 
W. Beenakker, W.L. van Neerven, R. Meng, G.A. Schuler, and J. Smith, NucL Phys. 
B351, 507 (1991). 

[22] A.D. Martin, W.J. Stirhng, R.S. Thorne, and G. Watt, Phys. Lett. B 652, 292 (2007) 
[arXiv:0706.0459j. 

[23] P.M. Nadolsky, H.-L. Lai, Q.-H. Cao, J. Huston, J. Pumphn, D. Stump, W.-K. Tung, and 
C.-P. Yuan, arXiv:0802.0007 [hep-ph]. 

[24] M. Czakon, A. Mitov, and S. Moch, Phys. Lett. B 651, 147 (2007), arXiv:0705. 1975 [hep- 
ph]; NucL Phys. B798, 210 (2008). arXiv:0707.4139. [hep-ph]: M. Czakon. ,arXiv:0803. 14001 
[hep-ph] . 



14 



